Dirt Softens Soap: Anomalous Elasticity of Disordered Smectics 
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We show that a smectic in a disordered medium (e.g., aerogel) exhibits anomalous elasticity, 
with the compression modulus B(k) vanishing and the bend modulus K(k) diverging as k — > 
0. In addition, the effective disorder develops long ranged correlations. These divergences are 
much stronger than those driven by thermal fluctuations in pure smectics, and are controlled by 
a zero temperature glassy fixed point, which we study in an e = 5 — d expansion. We discuss the 
experimental implications of these theoretical predictions. 
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The effects of quenched disorder on the properties of 
condensed matter systems continues to be a fascinating 
area of active research, which includes the study of dis- 
ordered superconductors fijl, charge density waves [2 
Josephson junction arrays and Helium in aerogel [4 
to name a few. Some of this attention has focussed [|](3|] 
on liquid crystals in the random environment of an aero- 
gel. While a complete picture of aerogel-confined liq- 
uid crystals is still being developed B], in this Letter 
we show that a smectic phase of these systems possesses 
strong anomalous elasticity, when subjected to an arbi- 
trarily weak amount of quenched disorder. This result 
has important experimental consequences. 

The anomalous elasticity of pure smectics was pre- 
dicted sometime ago [JtJ, and is characterized by bulk 
compressional and tilt moduli B(k) and K(k) which, re- 
spectively, vanish and diverge at long wavelengths (k — > 
0). This is a general property of all one-dimensional crys- 
tals in which the direction of the Id ordering wavevec- 
tor is chosen spontaneously. As a consequence of this 
spontaneous breaking of rotational symmetry, (a prop- 
erty of smectics but not of charge density waves) in such 
systems, compression can be relieved by smoothing out 
fluctuations (wrinkles in the smectic layers), so the effec- 
tive layer compressional modulus B(k) vanishes at long 
wavelengths. Similarly, in the presence of fluctuations, 
a bending of smectic layers necessarily leads to a com- 
pression, which implies that the effective tilt modulus 
if(k) diverges at long wavelengths. Unfortunately, this 
thermally driven behavior in pure smectics is difficult to 
observe experimentally, because the effect is very weak 
(logarithmic) in 3d. 

The main ingredients necessary for anomalous elastic- 
ity, namely spontaneously broken rotational invariance, 
and fluctuations, both exist even at zero temperature in 
quenched disordered smectics. In this Letter we demon- 
strate the existence of anomalous elasticity in quenched- 



disordered smectics which is significantly stronger than 
the marginal anomalous elasticity of thermal smectics, 
and is described by a zero-temperature fixed point that 
is perturbatively accessible in d = 5 — e dimensions. The 
elastic anomaly is much stronger in quenched disordered 
smectics because layer fluctuations are much larger, even 
at T = 0, than in a pure smectic at T > 0. 

One experimental signature of these divergences is in 
the smectic correlation length for smectics in aerogel, 
which has a different, universal dependence on the bare 
smectic elastic constants and the aerogel density than 
predicted by harmonic theory. Detailed predictions for 
these lengths can be found at the end of this Letter. 

Our model of a quenched disordered smectic starts 
with de Gennes' theory || for the smectic density field ip 
and the nematic director n, and includes a disorder field 
that couples to the nematic director, via 



d d r(g(i 



(1) 



where g(r) is a quenched random vector along the locally 
preferred nematic alignment. Such preferences can arise 
because, e.g., the nematogens may align with the local, 
randomly oriented aerogel strands. Below the nematic- 
smectic-A transition ip can be written as l^e 190 "^, with 
«(r) describing the local displacement of the smectic lay- 
ers from perfect periodic order. Furthermore, far below 
the nematic transition (i.e. inside the smectic phase) we 
can take n ks z + Sri, where z is the mean normal to the 
smectic layers, and _L denotes directions orthogonal to z. 
Integrating out the "massive" and <5n fields, which 
has the effect of replacing <5n — > V^u (i.e., the Higg's 
mechanism), and keeping only the most relevant terms, 
we obtain the elastic Hamiltonian that defines our model, 



H = 



K 



(Vi«) 



j(d z u-±(V^u) 2 f + h(r)-V ± u 



(2) 
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where h(r) = g z (r)g(r) is quenched random tilt disorder 
that for simplicity we take to be Gaussian, zero-mean, 
and completely characterized by 



hi^hjir') =A5 d (r-r')Si 



(3) 



The use of short-ranged correlations, even though the 
density of the (fractal) aerogel has long-ranged corre- 
lations, is justified, since the orientations of its con- 
stituent silica strands, beyond their microscopic persis- 
tance length, are certainly short-range correlated H|). 
We will focus here on the bahavior of the smectic in 
the (large) window of length scales between the intrin- 
sic orientational correlation length of the aerogel and the 
distance between disorder-induced smectic dislocations. 

The anharmonic terms included in Eq.^| are required by 
the underlying global rotational invariance of the smectic 
phase 0| , hidden by the spontaneous choice of the layers 
to stack along z-direction. 

To compute self-averaging quantities, e.g., the disorder 
averaged free energy, it is convenient (but not necessary) 
to employ the replica "trick" || that relies on the iden- 
tity log Z = lim„^o . This allows us to work with 
a translationally invariant field theory at the expense of 
introducing n replica fields (with the n — > limit to be 
taken at the end of the calculation). After replicating 
and integrating over the disorder h(r) utilizing Eq.^|, we 
obtain: 



1 r - 

H [u a ] = 2 / E 



K(V 2 ± u a ) 2 + B(d z u a - -(Vj.Ua) 



2T 



(4) 



a, /3=1 



from which the noninteracting propagator G a p(q) 
V^ 1 (u a (q)up(— q)) can be easily obtained, 

G al3 (q) = TG(q)S a(3 + Aq 2 ± G(q) 2 , 



(5) 



with G(q) = 1/ '(Kq^+Bql). The fluctuations associated 
with the disorder (the term in Eq|| proportional to A) 
are much larger, as q — > 0, than those associated with 
thermal fluctuations (the term proportional to T). 

We first attempt to assess the effects of the anhar- 
monicities, disorder and thermal fluctuations by perform- 
ing a simple perturbation expansion in the nonlinearities 
of H[u a ]. The lowest order correction SB to the bare 
elastic compressional modulus B comes from a part of 
the diagram in Fig.0. 




FIG. 1. Feynman graph that renormalizes the elastic mod- 
uli K, B and the disorder variance A. 



A simple analysis gives 
B 2 <■> 



SB 



f [TG(q) 2 + 2AqlG(q) 3 ] q 

Jq 
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dq z f> d q±_ 



_ x 2vr J (2tt)««-i (Kqi + Bq^f 



3 C d . 



16 5 - d \K 5 



1/2 



-5-d 



(6) 
(7) 
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where the constant C d = 2n d / 2 / ((27r) d r(d/2)) , we have 
dropped the thermal contribution that is subdominant 
to the disorder part, shown only the dominant contribu- 
tion for d < 5, and introduced a long wavelength cut- 
off L (defined by restricting the wavevector integral to 
q± > 1/L). The divergence of this correction as L — > oo 
signals the breakdown of conventional harmonic elastic 
theory on length scales longer than , which we define 
as the value of L at which \SB(^ L )\ = B. This definition 
gives: 



-,/2 



16(5 - d)K 
iC^B^A 



1/(5-0 



(9) 



which for physical 3d smectics is given by £j^ L = 

1 /2 

can also obtain a non-linear 
crossover length in the z-direction £J^ L = (£]v,l) 2 /A ~ 
K 2 1 'A, where A = (K/B) 1 / 2 , by imposing the infra-red 
cutoff in the z-direction. 

To understand the physics beyond this crossover scale, 
i.e. to make sense of the apparent infra-red divergences 
found in Eq.||, we employ the standard momentum shell 
renormalization group (RG) transformation. We sepa- 
rate the displacement field into high and low wavevector 
components: u a (r) — ii^(r) + w>(r), where has sup- 
port in the wavevector range Ae~^ < q± < A and A is an 
ultraviolet cutoff of order l/£jvL> integrate out the high 
wavevector part u^(r), and rescale the lengths and long 
wavelength part of the fields with r± — r' ± e e , z = z'e ul 
and u<(r) = e x u' a (r'), so as to restore the uv cutoff back 
to A. The underlying rotational invariance insures that 
the graphical corrections preserve the rotationally invari- 
ant operator {d z u a — i(Vj_ii Q ) 2 ), renormalizing it as a 
whole. It is therefore convenient (but not necessary) to 
choose the dimensional rescaling that also preserves this 
operator; the appropriate choice is x = 2 — lu. This rescal- 
ing then leads to the zeroth order RG flows of the effective 
couplings K{£) = ATe^" 1 " 1 ^, B{1) = _5 e (d+3-3u,y and 
(A/T)(£) = {A/T)e ( - d+1 -^ e . From these dimensional 
couplings one can construct two dimensionless couplings 
g x = (B/K 3 ) 1 / 2 and g 2 = A(B/AT 5 ) 1 / 2 , whose flow is 
given by g~i(£) = g\eC i ~ dS)l and g 2 (£) = f^e^ 5- ^. g~i is 
just the coupling that becomes relevant in d < 3 and was 
discovered in Ref. j?J to lead to anomalous elasticity in 
pure smectics. It is, however, only marginally irrelevant 



2 



in d = 3 and therefore only leads to a weak anomalous 
elasticity in physical 3d smectics. In contrast, the upper 
critical dimension d uc below which g 2 becomes relevant is 
due = 5, and leads to much stronger anomalous elasticity, 
that should be experimentally observable in disordered 
3d smectics. These observations imply that temperature 
is a strongly irrelevant variable near the disorder dom- 
inated fixed point. We will therefore set T = in all 
subsequent calculations. 

The integration over the high wavevector components 
of u a can only be accomplished perturbatively in non- 
linearities of H[u]. This perturbation theory can be rep- 
resented graphically; the graph giving the leading order 
corrections SB, 8K and 5 A (with the part diagonal in the 
replica indices a, (3 renormalizing K and -B, and the part 
independent of a, f3 correcting A) is shown in Fig.[j]. Jn| 
Evaluating it, and performing the rescalings described 
above, we obtain the following RG flow equations 



-^ = (d + 3-3»-- 92 )B 

dK ^ (A 1 _L 1 \K 



(10) 
(11) 



^™^ + l-^, 2 )(A/T), d2) 

where we have defined a dimensionless coupling constant 
g 2 = A( J B/AT 5 ) 1 / 2 C (i _iA <i - 5 , which obeys 



d92(i) _ 5 2 



(13) 



with e = 5 — d. As required the flow of g 2 is independent 
of the arbitrary choice of the anisotropy rescaling expo- 
nent u>. The RG flow Eq.|l3| shows that the Gaussian 
g 2 = fixed point becomes unstable for d < 5, and the 
low temperature phase is controlled by a stable, nontriv- 
ial, glassy T = fixed point at g 2 = 32e/5. 

The existence of this nontrivial fixed point leads to 
the anomalous elasticity, which we can calculate using 
the following matching approach. For this purpose it 
is convenient to use our RG results to evaluate the con- 
nected disordered averaged 2-point it(k) correlation func- 
tion G(k) oc (|u(k)| 2 ) - (u(k))(u(-k)). The power of 
the renormalization group is that it establishes a connec- 
tion between a correlation function at a small wavevector 
(which is impossible to calculate in perturbation theory 
due to the infra-red divergences) to the same correlation 
function at large wavevectors, which can be easily calcu- 
lated in a controlled perturbation theory. This relation 
for G(k) is 



G(k ± ,k z ,K,B,g 2 ) 



(14) 



e (3+«i-aO/ G ( kj y , kze ^ : K(£),B(i), g 2 (£)) , 

where the prefactor on the right hand side comes from 
the dimensional rescaling using the exact Ward identity 



X = 2 — uj, and we have traded in the disorder variable 
A for the dimensionless coupling g 2 . To establish the 
anomalous behavior of K, we look at k z = 0. We then 
choose the rescaling variable I* such that k±e e = A. 
We also choose k± sufficiently small such that gi{l*) has 
reached our nontrivial fixed point g 2 . Eliminating i* in 
favor of k i , we then obtain 



G(k ± ,0,K,B,g 2 ) 

A_ 
~k~ 



(15) 



3-\-d— uj 



G{k,Q,K{t),B{t),g* 2 ) , 



Since the right hand side is evaluated at the Brillouin 
zone boundary, it can be calculated perturbatively in the 
fixed point coupling g 2 . To lowest order we obtain 



G(k ± ,0,K,B,g 2 ) 



(A/fci) 



3+cZ— uj 



A 4 J fC(A/fc_ L ) 
1 

K(k ± )ki 



(d-l-o)+lff|/32) 



. (16) 



(17) 



where we integrated Eq.|ll| to obtain K(£*), and defined 
the anomalous tilt modulus which diverges at long length 
scales 



K(k ± ) = K(k ± /Ay 
with an anomalous exponent 



1 

32 
2 



1 



.92 



for d = 3 . 



(18) 

(19) 
(20) 



Similar calculations for the other coupling constants 
and other directions of k show that, in general, 



-t]K 



NL) 

fl(k) = B (k^ L ) m f B {KA z NL /{k^L) c ) , 
A(k) = A (k±i 



(21) 
(22) 



NL) 



fA(k z ^ N J(k^ NL ) c ) , (23) 



with the anisotropy exponent £ = 2 — (t)b + Vk)/^-, 
which would be 2 in the absence of anharmonic ef- 
fects, r\B = 3^1/16 = 6e/5 = 12/5 in d = 3, and 
VA = .92/64 = e/10 = 1/5 in d = 3. 

Of course, we do not completely trust the extrapola- 
tion of these small e results down to e = 2 [d = 3). 
However, since by definition dg 2 /d£ — at the nontriv- 
ial fixed point, this condition implies an exact relation 
between the anomalous exponents, 



5 - d + J] a 



2 



(24) 



which is obviously satisfied by the anomalous exponents, 
computed here to first order in e. This Ward identity 
between the anomalous exponents can be equally easily 
obtained from a self-consistent integral equation for the 
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u — u correlations functions, using renormalized wavevec- 
tor dependent elastic moduli and disorder variance, fl 

At length scales beyond £jj L and the elasticity 

and fluctuations of the disordered smectic are controlled 
by our new glassy fixed point. One of the important con- 
sequences can be seen in the layer fluctuations that can be 
observed in X-ray scattering experiments. For instance, 
layer displacement fluctuations along z arc described by 

C{z) = ((u(0 ± ,z)-u(0 ± ,0)) 2 ) , 

= f d d k 2(1 - cos(fc,z))A(k)fcj 

One can then naturally define the X-ray translational 
correlation length £^ , as the length along z at which 
C(z = = a 2 , where a is the smectic layer spacing. A 
simple calculation, using Eqs.|2~l|-p3| leads in 3d to 

g = {a/\f h K 2 /A = {a/X) 2h & L , (26) 

where 7 = (r/B +Vk)/C- Note that this X-ray correlation 
length is finite even as T — » 0. This result is consis- 
tent with the experimental observation || that the X- 
ray correlation length for smectics in aerogel saturates 
at some finite value at low temperatures. Note also that 
this length should be different for different smectics in 
the same aerogel, since B,K, and A will change from 
smectic to smectic. Since we expect A oc pa, the aerogel 
density the aerogel density dependence of Q could 
test the prediction of Eq.^. Likewise, the temperature 
dependence of £^ could be used to determine 7, since 
the bulk K(T) and B(T) that implicitly appear in Eq.|26| 
have temperature dependence that can be extracted from 
measurements on bulk materials. 

Note also that this correlation length is longer than the 
non-linear crossover length for A < a (i.e., for large B). 
For A >> a (small B), C(z) reaches a 2 before z reaches 
Qjl ' an d hence anharmonic effects are unimportant. In 
this case, the correlation length ££ can be determined 
in a harmonic theory (which amounts to evaluating the 
integral in Eq.^B] with K(k), B(k), and A(k) replaced 
by their constant (bare) values K, B, and A). This 
gives = a 2 BK/A = (a/A) 2 £^ £ , which is, reassur- 
ingly, much less than £,% L in the limit a « A in which 
we have asserted it applies. 

Although our entire discussion so far has focussed only 
on the effects of orientational disorder, we have shown |6J] 
that translational disorder (i.e., random pinning of the 
positions of the layers, is less important, at long wave- 
lengths (in d < 5), than the orientational disorder. Thus, 
the results described herein are directly applicable to real 
smectics, where both kinds of disorder are present. 

We have also ignored dislocations so far. Further de- 
tailed analysis || shows that in a harmonic theory dis- 
locations unbind for arbitrarily weak disorder. However, 



their effects are only felt on length scales greater than 
£,disi. = (a 2 /A)e*° A / A , where t is a parameter of or- 
der BK, which is much longer, in the weak disorder 
(A — > 0) limit, than the translational correlation length 
found here. Including anomalous elasticity, we find 
H a shorter ^disi. that nonetheless remains much longer 
than . Hence our predictions for and the anoma- 
lous elasticity studied here remain valid. 

In summary, we have found that the addition of disor- 
der to a smectic liquid crystal phase leads to anomalous 
elasticity, that is much stronger than in the pure sys- 
tem. This elasticity is controlled by a zero temperature 
glassy fixed point, perturbative in e = 5 — d. We expect 
that in contrast to marginally weak anomalous elastic- 
ity of pure smectics, the strong anomaly in the elasticity 
of disordered smectics, described here, should be readily 
observable. 

Note added: after this work was completely we learned 
of earlier work on the model Eq.|| in the context of 
spin glasses. However, this earlier work misses the es- 
sential physics of the nontrivial renormalization of disor- 
der, and therefore makes quantitatively and qualitatively 
wrong predictions of nearly all results; e.g. the anoma- 
lous exponents, disorder correlations, scaling exponent 
relations, etc., correct versions of which are given in our 
Eqs.|l2H24|. 
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